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I. INTRODUCTION 

In a scale larger than several Debye lengths, unmagnetized plasmas in the Universe have 
been assumed to be quasi-neutral in charge. This is because even a small charge imbalance 
would result in very large electric fields which would cause the very small time scale electro- 
static oscillation called plasma oscillation, and the damping or averaging of the oscillation 
over the mesoscopic time scale restores neutrality^*^. On the other hand, if we assume the 
ideal MHD condition on a magnetized, relativistically moving plasma, the charge density 
becomes significant, and this charge density is called the Goldreich- Julian density, per- It 
is noted the Goldreich- Julian density is required so that the electric field induced by the 
charge density vanishes in the comoving frame of the plasma. In such case, it seems that 
we cannot assume the charge quasi-neutrality. However, we point out here that the differ- 
ence between the net charge density of the magnetized plasma pe and the Goldreich- Julian 
density, Ape = pe — Pgj, plays the same role as the charge in the unmagnetized plasma. In 
this paper, we call the difference Ape the "free charge density", which induces the electric 
field observed by the comoving frame of the plasma. In the scale larger than several Debye 
lengths, the free charge density should tend to vanish because of the same reason for the 
charge quasi-neutrality of the unmagnetized plasma. This can be regarded as the gener- 
alization of the concept of charge quasi-neutrality to the magnetized, relativistic plasma. 
We call this concept "free charge quasi-neutrality". When the magnetic field is so strong 
and the thin plasma rotates so fast as assumed in a pulsar magnetosphere that the number 
density of the plasma particles is smaller than the Goldreich- Julian density divided by the 
elementary electric charge e, the free charge density becomes significant and the electric 
field observed by the comoving frame of the plasma remains. The parallel component of this 
electric field to the magnetic field accelerates plasma particles directly. In the region of the 
remaining one-direction electric field component along the magnetic field, plasma is swept 
by the electric field and the vacuum called "outer gap" appears^"^. Even in such vacuum, 
when plasma enters into it, the plasma is separated into positively and negatively charged 
fiuids and the two fiuids move to the opposite directions along the magnetic field to decrease 
the electric field component. Thus, in a strongly magnetized, relativistic plasma, the free 
charge tends to be canceled to keep the free charge neutrality. 

According to the above consideration, in a scale larger than several Debye lengths, it has 
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been assumed that the charged components of plasma move so that the electric field accel- 
erating the charged components decreases so as to restore the free charge quasi-neutrality. 
Here, we report a charge separation instability in an unmagnetized plasma rotating around 
a black hole, which will induce the free charge density and electric field exponentially. To 
investigate the charge separation of plasmas around the black holes, we use generalized 
GRMHD equations derived by Koide-. We present the linear analysis of the charge sep- 
aration of plasmas near Kerr black holes. We found the well-known plasma oscillation in 
the stable disk region outside of an innermost stable circular orbit (ISCO) around a black 
hole. On the other hand, in a circularly rotating plasma inside of the ISCO, we found an 
instability of the charge separation. The charge separation instability does not happen in 
the stable disk region outside of the ISCO. Furthermore, even in the unstable disk inside of 
the ISCO, the growth rate of the charge separation is smaller than that of the disk insta- 
bility. However, when the plasma density is much lower than the critical density (the very 
low plasma density case), the growth rates of the two instabilities become comparable and 
the charge separation instability becomes apparent. That is, due to the charge separation 
instability, the unstable disk falling into the black hole can be charged. 

In Section [Tll we present a linear analysis of charge separation in a stationarily rotating 
disk plasma around a Kerr black hole with a brief summary of the generalized GRMHD 
equations. The analysis shows an instability of charge separation in the plasma inside of the 
ISCO around a black hole. In Section IIIIl we summarize the results and discuss briefiy the 
astrophysical meanings of the charge separation instability. 

II. CHARGE SEPARATION IN PLASMA DISK AROUND BLACK HOLE 

We investigate the simplest process of charge separation in an unmagnetized disk rotating 
around a Kerr black hole. 

A. Brief summary of generalized GRMHD equations 

In this subsection, we briefiy summarize the generalized GRMHD equations; see Koide- 
in more detail. We investigate a charge separation using the generalized GRMHD equations 
of plasmas in the space-time, = {t, Ob ^ Ob ^ Ob ^) around a black hole where a line element ds 
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is given by ds^ = g^^dx^dx'' (Equations (18), (24), and (59) with Equations (25) and (58) 
of Koide-). Throughout this paper, except for a paragraph in Section [XTTt we use the unit 
system where light speed is unity and the energy densities of electric field E and magnetic 
field B are given by E"^ /2 and 5^/2 in the Minkowski space-time, respectively. 

First of all, we summarize the generalized GRMHD equations for a general case briefly 
as follows: 

V,{pW) = 0, (1) 
V^T^'^ = 0, (2) 

—V. i^Q^'^ = U^F^, -r]{J^- p'eUn + -^^^(A/ip - Ap) 
ne \ ne J 2ne 

^^-rF^^, + r]p'^eU^, (3) 



ne 

and Maxwell equations 

V, *F'"' = 0, (4) 
V^F^"" = J'', (5) 

where the energy-momentum tensor T^'^ and "charge-current density tensor" Q^'^ are given 
by 

[neY ne 

_ i^^M-F'^^F.A, (6) 



^ - pW ~ W 



' hX Ahi 

— (Wr + J^U") + 2Ah^Uf'U'' - ^^J^J 
ne (ne)'' 



(7) 



Equation ([3]) presents the general relativistic generalized Ohm's law. In Equation ([3]), the 
left-hand side expresses the inertia effect and transport of kinetic energy and momentum 
of the current, the first two terms of the right-hand side correspond to all terms of the 
"standard" Ohm's law with resistivity t], the third term represents the thermo-electromotive 
force, the forth term expresses the Hall effect, and the last term comes from the equipartition 
of the thermalized energy due to the friction force between the two fluids. Here, Pg = —J'^J^ 
is the charge density observed by the local rest frame of the plasma and G is the rate of 
equipartition with respect to the thermalized energy due to friction (for detail, see Appendix 
A of Koide ^). We follow the notations used by Koide- with respect to physical variables 
except that we use Q^^ instead of K^^ . Here, we used the two- fluid model, where we 
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assumed the plasma consists of positively charged particles with charge e and mass m+ 
and negatively charged particles with charge — e and mass m_ (Appendix [A]) . We used the 
typical mass of a plasma particle m = m+ + m_, normalized reduced mass fi = m^m^/m'^, 
and normalized mass difference A/i = (m+ — mJ)/m. The variables p, h\ p, n = p/m, Ap, 
and Ah"^ are the mass density, enthalpy density, pressure, number density, pressure difference 
of two fluids, and difference of two fluid enthalpy density. Furthermore, V^, If^, and J'^ 
are the covariant derivative, 4-velocity, and 4-current densities, respectively, and F^^^, is the 
electromagnetic strength tensor and *F^u is the dual tensor of Here, the electric field is 
given by Ei = FiQ and the magnetic field is Fij = J2k ^ijkBk {^ijk is the Levi-Civita symbol), 
where the alphabetic index (i, j, k) runs from 1 to 3. We also use the variables related to 
the enthalpy density, 

h} = h^~ ApAh\ Ah^ = Aph) - l^^A/zt. (8) 

It is noted that Equation ([5]) yields the equation of continuity with respect to the current, 

= 0. (9) 



We assume that off-diagonal spatial elements of the metric g^^iy vanish: gij = {i ^ j). 
Writing the non-zero components by g^Q = —h^, gu = hf, giQ = goi = —hfui, we have ds"^ = 
g^^dx^dx" = —hldt"^ + J2^=i [^f l*^^*)^ ~ 2h^ujidtdx^]. When we define the lapse function a 

1 /2 

and shift vector /3* by a = [/iq + X]i=i {hi^i) ] , /3* = the line element ds is written 
by ds^ = —a^dt^ + ~ afS^dt)"^. We also have g = —{ahih2h^Y . Using the "zero- 

angular-momentum observer (ZAMO) frame" x^, where the line element ds is given by 
ds^ = —dv + Ylii^^y — V/j.i^dx^dx'^ , we have the 3+1 formalism of the generalized GRMHD 
and the Maxwell equations. As for equations including only derivatives of contravariant 
vectors or anti-symmetric 2nd rank tensors A^"^, we obtain their 3+1 formalism easily 
using Equations V^A" = -^d{^/^A'') or V^A^'' = -^d^^/^A^""). With respect to any 
equation including a term of derivative of the symmetric 2nd rank tensor. 
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the 3+1 formalism is given by 



d 



d 



dt hih^h'i ^ dx^ 



ahih2hz 
hi 



\ / hj ox^ 



j,k j,k 



d_ 
di 



5* + 



d 



/ll/l2^3 dx^ 



ahih2h3 



+ 



1 da 
hi dx' 



S 



00 



(11) 



(12) 



where dj = -^ff and a.j = ^^(a/^O- 



B. Linecir analysis of charge separation in stationary disk 

For simplicity, we consider a plasma of a stationary thin disk rotating around a Kerr black 
hole with zero pressure {p = Ap = 0). The space-time x'^ = (t, r, 6, (j)) around the Kerr black 
hole with a mass M and rotation parameter a is given by the metrics, /iq = (1 — 2rgr/S)^/^, 
hi = V^/A, /i2 = VE, /is = V^/E sin 6*, uj^ = 2rlar/A, and = (i = 1,2). Here, 
rg = GM is the gravitational radius {G is the gravitational constant), A = r^ — 2rgr+(arg)^, 
E = + (arg)^cos^^, and A = {r^ + (arg)^}^ — A(arg)^ sin^ ^. In this metric, the lapse 
function is a = a/AS/A. The Schwarzschild radius of the black hole is given by = 2rg. 
The 3-velocity of the circularly rotating disk observed by the ZAMO frame, called the Kepler 
velocity, Vk, is given by the quadratic equation 

aG.iVi + (a/J^Gai + aaO^x + 7-^ = 0. (13) 

All (Jf 

In investigating linear behavior of charge separation in the stationary disk, we assume charge 
separation is weak, \p'^\ = l — J^J^l ^ en. Then, we can use an approximation of the enthalpy 
density and enthalpy difference density as 

h'^ ^ mn, (14) 

A/i^ ~ ^P'e^ (15) 

h^ ~ mn H Pg, (16) 

A/i" ^ mnAii + - — —mp'^, (17) 
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because of n± ^ n ± ^Pg- The generalized GRMHD equations reduce to 



V^ipUn = 0, (18) 



ner V en (ne 



J'F^, (19) 
(20) 



where Up = a/ (ne)^/ (pp) = ^Jne^ j (fim) is the plasma frequency, and Q^'^ is approximated 
by 

Q'^'^ ^ u^'r + - p;,f/^f/^ (21) 

To perform the linear analysis of the charge separation in the hydrostatic equilibrium plasma 
rotating around the Kerr black hole, we consider only the perturbation with respect to the 
static electric field, 

> = >, pe = Pe (22) 

Fio = Ei, Fij = 0, (23) 

where the tildes indicate the infinitesimally small variables, and we do not consider pertur- 
bation to the hydrostatic equilibrium, 

p = p, n = n, U^' = U^. (24) 

In the ZAMO frame, the 4-velocity is given by t/° = 7k = (1 - Vl)-^/"^, tj^ = U'^ = 0, 
= 'JkVk- The linear analysis requires the Ohm's law (Equation ([3])), the equation 
of continuity about current (Equation ([9])), and the Gauss law of electrostatics (temporal 
component of Equation ([5])), 

\V.Q^'' = U'^Fi: - r/[> - (p1 + p'^e)U% (25) 

Vj" = 0, (26) 
V^F^'' = J°. (27) 

Using Equations (I25l) - (p7j) and (fT2|) . we obtain the 3+1 formalism of the Ohm's law, equation 
of continuity about current, and Gauss law for the electric field (see also Equations (63) and 
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(67) of Koide-^), 



1 ^ d 



d ( ahxh2h-^ 



1 da ~ 



hj dx 



d 



d fhih2h3jj 



Pe 



^i^2^3 dx^ \ hn 
J 

. hih2h3dx^ V ^ 



(28) 
(29) 

(30) 



Here, we used the following approximation, 



g™^7K(2pe-7KPe), 

Q^o = Q0^^7KJ* + (pe-7KP;)f/\ 



(31) 
(32) 
(33) 



where Q'^^ can be regarded as modified charge density and Q'^ corresponds to the modified 
current density. When we use the relation 

= J^Uu = -fpe - J ■ U ^ 7kPo -tJ J, 

where J = {J^,J^,J^) and U = {U^,U^,U^) are the 3-current density and 3-velocity, 
respectively, we have 



g""~7K[(l-f/^)Pc + 7Kf/KJ1, 



(34) 
(35) 



Because in the present linear analysis, we can assume the quasi-charge neutrality, thus we 
have an approximation of p'^G as. 



^|j/|2 

2ne ' 



(36) 



where |J'p = (7 — l)pg — 27PcC/- J+ | + (t/- J)^.- Then, when J' is infinitesimally small, 
we have p'^Q ~ 0. 
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Here, we assume the perturbation is symmetric with respect to the polar axis and the 
equatorial plane, d/dx^ = 6^^d/dx^. Then, the equations with respect to the perturbation 
of the charge separation at the equatorial plane are as follows: 



2 



7k 

dpc 1 d 



|-Ji = [«G3i(2\/k + + asi]{VKPe - J') + ^(7K^i - VJ'), (37) 
at 7k 

^3 - V {ikJ' - CkP.)] , (39) 

To derive Equation fl57|) . we used Equation f[T^ . For simplicity, we assume the wave length 
of the perturbation is much smaller than the characteristic scale-length of the metrics and the 
Keplerian rotation around the black hole, ~ rg, and we put the perturbation is proportional 
to exp{ikr — iujt). Here, we note that we have to consider the relation between the derivatives 

^ ^ [ahhlUKJ') - , ] , ^{ahhlJ') = (42) 



aUKhih2h\dr ahih2h^dr Lk dr ' 

where Lk = h^U^ is the specific angular momentum of the Keplerian disk. Finally, we 
obtain the dispersion relation of the charge separation in the plasma disk rotating circularly 
around the Kerr black hole, 

where rj' = Upij/jii, 

9K = ^ [aGsii2VK + f) + asi] = -^^^ loga - ^("Ggi/?^ + ^31), (44) 
/ii7^Lk dr 

Here, it is noted that the stability condition of the accretion disk is given by dL]<^/dr > 
(Ak > 0), and the condition dL^/dr = (Ak = 0) yields the radial coordinate of the ISCO, 
r = nsco- 

In the case of zero resistivity (77 = 0), Equation ( l43l) yields the dispersion relation of the 
charge separation in the plasma disk as 

= a\ul - 2gKAK) ^ u;l (46) 
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In the region inside of the ISCO, when 



< v/2(7kAk, (47) 



the charge separation becomes unstable. However, it is noted that the disk is unstable 
with the growth rate 7disk = Q^^f^QK^Ki which is larger than the growth rate of the charge 
separation instability. In the very low plasma density case (wp <^ 7disk), the growth rate of 
the disk instability and the charge separation instability become comparable. Then, in this 
situation, the charge separation instability may appear in the disk falling into the black hole 
and may make the black hole charged. On the other hand, in the other usual plasma density 
case, the charge separation instability is forbidden or inhibited behind the disk instability. 

To investigate the effect of resistivity, we consider a solution for a very weak resistivity 
limit, arj' ^ Wq. In this approach, we treat the difference Aw = — is an infinitesimal 
variable which is comparable to {arj' /ujI)ujq. The dispersion relation ( H3|) yields 

2w^Aw + iri' [2ujI - a^wj] = 0, (48) 

and we have the solution 

oi^v' / 9 . N / 9 . M/9 «^V(^n ~ 45(kAk) , , 

It is noted that even in the case of > (stable state in the zero resistivity case), Up—AqkAk 
can become negative in the unstable region (Ak < 0), which means the resistivity induces 
the charge separation instability, while in the outside of the ISCO, r > risco; the instability 
is forbidden. On the other hand, in the case of Uq < 0, the resistivity stabilizes the charge 
separation instability as shown by the last term of the right-hand side of Equation fH^]) . 



III. DISCUSSION 

We have shown the charge separation instability of the circularly rotating plasma inside 
the ISCO (r < risco) around the Kerr black hole. This instability is forbidden in the stable 
disk outside of the ISCO. Furthermore, even in the unstable disk region inside of the ISCO, 
the growth rate of the charge separation instability is smaller than that of the disk instability. 
However, when the plasma density is much lower than the critical density, growth rates of 
the charge separation and disk instabilities become comparable and the charge separation 

10 



Charge separation instability around black hole 



instability becomes apparent. Then, the charge separation instability makes the disk plasma 
falling into the black hole charged and may eventually charge the hole. The critical plasma 
density is expressed as 

^crit = — ttQkJ^k, (50) 

where the critical density ncrit is given by Equation ( H71) with = ncrite^/(/im) in the zero 
resistivity case. With respect to the resistive case, Equation (H^ indicates that the weak 
resistivity makes the charge separation instability happen easier in the stable region of the 
charge separation instability inside of the ISCO, while in the unstable region of the ideal 
MHD situation, the resistivity stabilizes the instability. In a case with resistivity, the growth 
rate of the charge separation instability is also smaller than that of the disk instability inside 
of the ISCO. Then, unless the plasma density n is much less than the critical density ricrit 
so that the growth rate of the charge separation instability is comparable to that of the 
disk instability, the charge separation instability is hidden behind the disk instability inside 
of the ISCO and is also forbidden in the stable disk around the astrophysical black holes. 
However, we emphasize again that when the plasma density is much less than the critical 
density, n ^ ^crit, the charge separation instability may be apparent as the charged disk 
falling into the black hole. 

The mechanism of the charge separation instability is explained by the following schematic 
picture. Consider that the accretion disk is composed of two disks of positively charged 
particles (+ disk) and of negatively charged particles (- disk) (Fig. [I]). Both of the purely 
charged disks (± disks) are unstable in the region r < risco- However, even in the unstable 
region, when one disk (± disk) falls and the other disk (=p disk) shifts outward, the electric 
field is induced and tends to suppress the charge separation. When the electric field is 
strong enough, it causes plasma oscillation. On the other hand, when the electric field can 
not become so strong (in the case of n ^ ^T-crit)? the suppression is not strong (effective) 
enough and the charged disks separate each other increasingly so that the charge separation 
is induced exponentially. In the outer region of the ISCO, r > risco, the plasma oscillation 
is always induced by the charge separation because of the stability of the ± disks. Using this 
intuitive picture of the charge separation instability, we suggest that the charge separation 
instability with small wave number k is also possible in the unstable disk around the rotating 
black hole, when the plasma oscillation frequency is small enough and the plasma density 
is sufficiently small. This picture also suggests that in a magnetized plasma disk, the similar 
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instability of charge separation would be caused to break the free charge quasi-neutrality 
when the purely charged disk is gravitationally unstable and the plasma density is low 
enough. This charge separation grows exponentially until the nonlinear effects begin to 
suppress the instability. For example, the disk charged by the charge separation instability 
falls into the black hole and charges the black hole. The electric field of the charged black 
hole will suppress the fall of the charged disk and induce the fall of the oppositely charged 
plasma around the black hole. 

The charge separation instability may also happen in the pulsar magnetosphere when 
the radius of the ISCO is larger than the central star radius. However, when the magnetic 
field is extremely strong, the free charge due to the charge separation instability would 
not be induced because the magnetic field suppresses the gravitational instability of the 
disk. Plasma dynamics of the pulsar magnetosphere outside of the relativistic star were 
investigated with the Schwarzschild metric by Henriksen and Rayburn^. They discussed the 
net charge separation and the instability, while they did not perform the linear analysis of 
the charge separation in bulk plasma. They considered the large-scale net charge separation 
within the free charge quasi-neutrality and the charge separation instability produced by the 
current driven plasma turbulence which leads charge fluctuation on small scales less than 
the Debye length. 

Here, we estimate the critical density of the charge separation instability (Equation ( 150|) ) 
for an individual astrophysical object. For a rough estimation of the instability condition, 
we consider the cases of the Schwarzschild black holes with the mass Mbh = 0). The 
critical density of the instability is expressed in the MKSA system of units (SI unit) by 

_ fimrs 3rs - r 

where /iq is the magnetic permeability in vacuum. We estimate it at r = 2rs to get 

<« = ^W^ = 2-5xl""('-=H)-^^^ h-'] = 2.8. 10^ (^yt^, (52) 

me 2''^fiQe^r^ rrie \ Mq / rrie 

where rrte is the electron mass. If we assume the electron-proton and electron-positron 
plasmas, we have fim/rrie = 1 and fim/me = 1/2, respectively. In this paragraph, we set 
fim/nie = 1. We use values of the typical density n of accretion disks around black holes 
of individual objects listed up by Koide-. When we consider the active galactic nucleus 
(AGN) of M87 whose central black hole mass is Mbh = 3 x lO^M© (rs = 9 x W'^m^, the 
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critical density of the charge separation instability is estimated as n*^,;^ = 3 x 10~^^ [m~^]. 
This value is extremely small compared not only to the estimated value at the accretion disk 
around the black hole of = 8 x 10^^ but also to the averaged particle number density 

in the extragalactic region, riu ~ 10 [m^^].^*^ In the case of Sgr A*, a supermassive black 
hole in Galaxy, whose central black hole mass is Mbh = 4.4 x lO^M© (rs = 1.3 x lO^^m)^, 
the critical density is n*^,;^ = 1.4 x 10~^[m~^]. This is also extremely small compared to 
the value at the accretion disk, n = 9 x 10^^ [m^'^]. As we estimated above, the charge 
separation instability hardly happens around the supermassive black holes. In the case of 
black hole X-ray binaries, for example, the micro-quasar, GRS191 5+105, whose central black 
hole mass is Mbh = 14M0 (rs = 4.2 x lO'^m)-'^^, the critical density is n^j.-^^^ = 1.4 x 10^ [m^^]. 
This density is much smaller than the density at the accretion disk around the black hole, 
n = 4 X 10^^ [m~^]. However, a density smaller than the critical density may be realized in 
a low density disk around a single stellar-mass black hole, for example. In the case of a very 
low density disk, the charge separation instability may be caused, while it would be difficult 
to observe because of its low activity of the thin disk. 

The charge instability causes the falling of the charged plasma into the black hole and will 
charge the black hole. Here, we estimate the influence of the electric fleld of the charged black 
hole to the ambient plasma in the following case. We assume that the one-component charged 
fluid with the height ~ rs, the inner radius runner ~ 2rs, the outer radius router ~ 2rs + L, 
and the density n ~ n*^.^^ falls due to the charge separation instability and is swallowed by 
the black hole. The charge contained in the one-component fluid is Q ~ 87rrgL?2*j.j^e. The 
electric and gravitational forces which act on a charged particle with the mass m_ and the 
charge — e located at r = r are ~ Fgrav ~ <j "^^2°" , respectively. We estimate 

the disk scale L where the electric and gravitational forces become comparable: ~ -Fgrav 

2 

It yields L ~ /""V^ . Using Equation (!52l) . we have the simple expression, 

L ^ 2^/2^ ~ 2.8, (53) 
rs m+ 

where the ratio does not depend on the black hole mass, Mbh- This suggests that the small 
disk with the scale of the Schwarzschild radius supplies so much charge to the black hole 
through the charge separation instability within short time scale of the instability that the 
charged black hole influences the plasma dynamics around it. This electric fleld will suppress 
the further charging of the black hole. 
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It is true that the charge separation instability is caused only in the unstable disk region. 
However, the current induced by the instability may reach the stable disk region because 
of the inertia of the current. The current supplies the net charge in the stable disk around 
the black hole. Thus, the charge separation due to the instability may induce a distinctive 
drastic phenomena in the plasma of the stable disk regions around black holes, while the 
net charge may be canceled by the charge supply from the outer disk. For example, the 
gravitational magnetic reconnection can be induced by the charge separation^. If the charge 
separation is significant, it also causes a strong electric field. This strong electric field may 
accelerate particles which can explain high energy cosmic rays as well as in the "outer gap" 
of the pulsar magnetospher^>5,. The cause of these distinctive phenomena of plasma around 
black holes will be clarified with more detailed analysis of the generalized GRMHD equations 
and numerical simulations. The numerical simulations of the generalized GRMHD beyond 
the ideal GRMHD^ii^ and the resistive relativistic MHD with the acausal relativistic Ohm's 
law^ would provide a useful and essential tool for such analysis. This is our next subject. 

In the last paragraph of this paper, we suggest a radiation process associated with the 
charge separation instability. The charge separation instability does not depend on the wave 
number k in the zero resistivity case as shown by Equation ( l43l) . Then, very strong electric 
field with a very large wave number can be caused through the instability. When a high 
energy charged particle passes through the strong electric field, the particle is accelerated 
and decelerated reciprocally, and emits strong radiation. This radiation should be detected 
around the inner edge of disks of black holes. 
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FIG. 1. Schematic picture of the charge separation instability of the plasma in an unstable region 
around a black hole (r < nsco)- The neutral disk consists of positively /negatively charged disks 
(+/- disks); when the electric field due to the charge separation is not strong enough to suppress 
the charge separation, for example, the positive disk shifts outward and the negative disk falls into 
the black hole exponentially in the unstable disk region. The reverse is also true. This shift and 
falling cause the exponential charge separation. 



Appendix A: Basis of generalized GRMHD equations 



The generalized GRMHD equations are derived from the general relativistic two-fluid 
equations of plasma, which is composed of positively charged particles with charge e and 
mass m+ and negatively charged particles with charge e and mass m_-. The variables of 
the generalized GRMHD equations are deflned by the average and difference of the variables 
of the two-fluid equations. We note a variable of fluid composed by the positively charged 
particles with a subscript "+" and that of negative fluid with "-" . We write the variables of 
two fluids as: n± is the number density, 7^ is the Lorentz factor observed by the center of 
mass frame of the charged fluids, p± is pressure, h± is the enthalpy density, U'^ is 4-velocity. 
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The variables of the generalized GRMHD equations are based as follows: 
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n = —, 
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P = P++P-, 
Ap = p+ 
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We also use the variables with respect to the enthalpy density: 
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Incidentally, we sometimes assume the plasma consists of two perfect fluids with the equal 
specific heat ratio, F. The equations of states are 

r 



n 



Ah'^ = 2/xmn^ 



n„ 2(r — 1) 1^ 
1 1 

1 1 \ / 1 



n, n_ 



(All) 



2(r — 1) \ \m+n\ m-n_ 



P 



+ 



Ap 



where 



1/2 



(A12) 



(A13) 



em \em. 

corresponds to the particle number density of each charged fluid (see Equations (74) - (78) 
of Koide-^l]. 



^ Equations (76) of Koide^ contain several typographical errors. Here, we correct them in Equation (jAlSp . 
The variable p± in Equation (76) of Koide— is given by p± = mn±. 
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The energy-momentum tensors Ti^" and the charge-current density tensor Qi^" are given 

by 

Qi^- = ^^( J-T!r - —rr] , (A15) 
AiT \m+ m_ / 

where T^'^ — g^'^p± + h±U^U^ are the energy-momentum tensor of the two fluids and 
T^^ = Fi'^F"" - Igi'^F^^F^x is the Maxwell stress tensor. 

The relativistic two-fluid equations come from the continuity equations of particle number 

and conservation law of energy and momentum: 

V.(n±C/^) = 0, (A16) 
V,{h±U^U^) = -V'pi ± en±U^^ ± R^, (A17) 

where i?'* is the frictional 4-force density between the two fluids. It connected with the 
resistivity and current as 

R^" ^ -rjneiJ" - p'^{l + ep^"]. (A18) 
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